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Abstract 

We study the light-element abundances in super symmetric model where the right- 
handed sneutrino is the lightest superparticle (LSP), assuming that the neutrino masses 
are purely Dirac-type. In such a scenario, the lightest superparticle in the minimal 
supersymmetric standard model sector (which we call MSSM-LSP) becomes long-lived, 
and thermal relic MSSM-LSP may decay after the big-bang nucleosynthesis starts. 
We calculate the light-element abundances including non-standard nuclear reactions 
induced by the MSSM-LSP decay, and derive constraints on the scenario of right- 
handed sneutrino LSP. 



With the precise astrophysical observations, it is now widely beheved that about 23% of 
the energy density of the present universe is due to dark matter (DM) The existence 
of dark matter, however, raises a serious question to particle physics because there is no 
viable candidate for dark matter in the particle content of the standard model. To solve this 
problem, many dark-matter models have been proposed so far. 

In constructing dark-matter model, it is important to understand how dark matter was 
produced in the early universe. In many cases, the thermal freeze-out mechanism is adopted 
to produce dark matter particle in the early universe; then, dark matter particle, which is 
in thermal bath when the cosmic temperature is higher than its mass, freezes out from the 
thermal bath when the cosmic temperature becomes low. 

However, the freeze-out scenario is not the only possibility to produce dark matter particle 
in the early universe. Even if the dark matter particle is very weakly interacting so that it 
is never thermalized, it can be produced by the decay and scattering of particles in thermal 
bath. In particular, if the interaction of dark matter is dominated by renormalizable ones, 
dark-matter production is most effective when the temperature is comparable to the mass 
of parent particle which produces dark matter via the decay or scattering. Thus, if the 
reheating temperature after inflation is higher than the mass of parent particle, the relic 
density of the dark matter becomes insensitive to the cosmic evolution in the early stage. 

Such a scenario was originally proposed in [2], where the right-handed sneutrino in 
supersymmetric model is shown to be a viable candidate for dark matter. In [2], it was also 
shown that, if z/jj-DM is dominantly produced from the decay and scattering of superparticles 
in thermal bath, the primordial abundance of is determined when the cosmic temperature 
is comparable to the masses of superparticles. Then, recently, more general discussion of 
such a scenario has been given in [3J, where a variety of candidates for such very weakly- 
interacting dark-matter particles have been also considered. 

If a very weakly interacting particle is dark matter, it is often the case that a long-lived 
particle (with lifetime longer than ~ 1 sec) may show up, which may spoil the success of the 
standard big-bang nucleosynthesis [H El |2l [6], [3] . This is indeed the case where the right- 
handed sneutrino is the lightest superparticle (LSP) and is dark matter. If a right-handed 
sneutrino is the LSP, the lightest superparticle in the minimal supersymmetric standard 
model (MSSM) sector (which we call MSSM-LSP) decays into (and i?-even particles) via 
very small neutrino Yukawa interaction. Then, decay of the MSSM-LSP after the big-bang 
nucleosynthesis (BBN) epoch may affect the light-element abundances. Thus, it is important 
to check the BBN constraints on the scenario. 

In this letter, we consider the case where a right-handed sneutrino is the LSP, assuming 



Dirac-type neutrino masses We study the light-element abundances in such a case 

in detail, and derive BBN constraints on the mass and lifetime of the MSSM-LSP. We 
also comment on the implication of the sneutrino LSP scenario on the ^Li overproduction 
problem. 

First, we discuss the model framework that we consider in this letter. The superpotential 



#ipQ]- related topics, see also [7]. 



1 



is written as 



W = WussM + VvLHupR, (1) 

where PFmssm is the superpotential of the MSSM, L = {vl^^^l) and Hu = {H^,H'^) are 
left-handed lepton doublet and up- type Higgs doublet, respectively. (In this letter, "hat" is 
used for superfields, while "tilde" is for superpartners.) Generation indices are omitted for 
simplicity. In this model, neutrinos acquire their masses only through Yukawa interactions 
as nil, = yu{H'^) = Z/i^'i^ sin /?, where v ~ 174 GeV is the vacuum expectation value (VEV) of 
the standard model Higgs field and tan jd = {H^) / {H^)- Thus, the neutrino Yukawa coupling 
is determined by the neutrino mass as 

l/^sin/? = 3.0 X 10"^^ X ^- ^ . (2) 

^ ^ V2.8 X 10-3 eVV 

Mass squared differences among neutrinos have already been determined accurately by neu- 
trino oscillation experiments. In particular, the K2K experiment suggests [Am^jatom — 
(1.9 — 3.5) X 10"^ eV^ [8]. In the following discussion, we assume that the spectrum of 
neutrino masses is hierarchical, hence the largest neutrino Yukawa coupling is of the order 
of 10"^^ unless otherwise mentioned; we use i/i, = 3.0 x 10-^^ for our numerical study. (We 
neglect effects of smaller Yukawa coupling constants.) For our study, it is also necessary to 
introduce soft supersymmetry (SUSY) breaking terms. Soft SUSY breaking terms relevant 
to our analysis are 

Aoft = 'lim^BB + m^WW + h.c.) - Ml^L - Mj^u^Vn + {A,LHuU% + h.c), (3) 

where B and W are Bino and Wino, respectively. We parameterize A^, by using the dimen- 
sionless constant 

A^ = a^y^Mi. (4) 

Notice that is a free parameter and, in gravity-mediated SUSY breaking scenario, for 
example, is expected to be 0(1). The A^^-term induces the left-right mixing in the 
sneutrino mass matrix, through which the MSSM-LSP decays in the present case. In the 
calculation of mass eigenvalues, however, the mixing is negligible because of the smallness 
of neutrino Yukawa coupling constants, and we obtain 

m?,c.M| + i,4coB2ft (5) 

where mz is the Z boson mass. Here and hereafter, we assume that all the right-handed 
sneutrinos are degenerate in mass for simplicity. In the numerical study, we take the following 
model parameters: m^j^ = 100 GeV, tan/? = 30, and rrih = 115 GeV (with nih being the 
lightest Higgs boson mass). In addition, the Wino mass is related to the Bino mass using 
the GUT relation. 
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In the early universe, right-handed sneutrino is never thermahzed because of the weakness 
of neutrino Yukawa interaction. Although it is decoupled from thermal bath, right-handed 
sneutrino can be produced in various processes; (i) decay or scattering of MSSM particles 
in thermal bath, (ii) decay of MSSM-LSP after freeze-out, and (iii) production in very early 
universe via the decay of exotic particles (like gravitino or inflaton). Thereafter, we donate 
the contribution of each process as, j^f.™'"'™'''' ^ and (^(."^""-^s^^''' order. Primarily, 

right-handed sneutrino is produced through neutrino Yukawa interaction (and the left-right 
mixing of sneutrino) dominantly in the following decay processes: —>■ vrv., — »• vrI^ , 
vl — ^ vrU^ vl t^rZ , II — ^ i^RW~ , B — >• pru, — >• vrv, and — * urI^ . In the previous 
work, it was shown that right-handed sneutrino can be adequately produced to become dark 
matter when the masses of left- and right-handed sneutrino are degenerate at 10 — 20 % with 

< 3, or in a case of larger a^, without degeneracy [2]. It is also mentioned that enhancement 
of right-handed sneutrino production is possible with larger neutrino Yukawa coupling if we 
consider the case where neutrino masses are degenerate. Giving an eye on the thermal bath 
again, the MSSM-LSP decouples from thermal bath and its number freezes out in the same 
manner with usual MSSM, while the number of the other MSSM particles is suppressed by 
Boltzmann factor in this epoch. However, relic MSSM-LSP, which is assumed to be the 
next-to-the-lightest superparticle (NLSP) in this letter, decays to right-handed sneutrino 
through neutrino Yukawa coupling in the late time. In this process, the contribution to the 
abundance is given as 

(^(F.o.) _ rrij^E^ n(F-O-) 

^^UR ~ ^ "nlsp' y^i 

where mNLSP is the mass of the NLSP and ^^nlsp would-be density parameter of 

the relic NLSP (for the case where it does not decay into vr). Lastly, we mention that 
there might be a possibility that right-handed sneutrino is produced directly from an exotic 
particle in the very early universe. The abundance of the expected right-handed sneutrino 
is model-dependent, and we do not discuss further detail of specific model. In this letter, we 
consider the scenario that right-handed sneutrino produced in these processes becomes dark 
matter. We do not specify which is dominant process to produce right-handed sneutrino. 

If vr is the LSP, it is always the case that the MSSM-LSP becomes long-lived. Because 
some amount of relic NLSP always exists in the early universe, they may cause serious 
problem in BEN; if the relic MSSM-LSP decays during or after the BEN epoch, energetic 
charged and/or colored particles are emitted; they cause the photo- and hadro-dissociation 
processes of light elements, which may spoil the success of the standard EEN scenario. In 
the following, we consider three typical candidates for the MSSM-LSP; Eino 5, left-handed 
sneutrino z/^, and lighter stau f, and study how the z/^-DM scenario is constrained by the 
EEN. 

In the Eino-NLSP case, the Eino dominantly decays as B ^ vrv (and its CP-conjugated 
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process) and its decay rate is given 
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where gi is the U{1)y gauge couphng constant and, for the process x z/^y, /3f is given by 

1 



(it 



mi 



r 4 o/ 2 I 2\ 2 I / 2 2\2i 



with and being the masses of the particles x and respectively. When or f is 
the NLSP, the NLSP decays by emitting weak- or Higgs-boson if kinematically allowed. The 
decay rates for those processes are given by 
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where mf is the stau mass and Of is the left-right mixing angle of stau. (The lighter stau is 
given by f = TRCOs^f -f- fLsin6'f.) If the two-body processes are kinematically blocked, the 
slepton-NLSP decays into three-body final state as vl ^rII and f — *• I'nf f (with / and 
/' being standard- model fermions). 

Now, we are at the position to discuss the BBN constraints on zz/j-DM scenario. We start 
with the case where Bino is NLSP. As we have mentioned, the Bino-NLSP dominantly decays 
as B ^ vrv. Since vr and v are (very) weakly interacting particles, the BBN constraints 
are not so severe if this is the only possible decay mode. However, B may also decay as 
B uruZ'-*^ and i>RlW^*\ where Z^*) and W^*^ are on-shell or off-shell Z and W bosons 
(where the "star" is for off-shell particle), respectively, while / is charged lepton. Then, 
through the decay of Z^*^ and W^*^ , quarks and charged leptons are produced. Even though 
the branching ratio for such processes are phase-space suppressed, they produce sizable 
amount of hadrons which may significantly affect the light-element abundances. Thus, in our 
analysis, effects of those decay modes are taken into account in deriving the BBN constraints. 
The light-element abundances also depend on the primordial abundance of the NLSP, and 
we adopt the abundance of Bino in the focus-point (or co-annihilation) region 



y(focus) ^ ^ 13 f 

s VlOOGeV/ 



(12) 



"^^In this letter, we consider the case where the Gaugino-Higgsino mixing is so small that its effect is 
negligible. 

^•^If the Bino is the NLSP, its primordial abundance strongly depends on the MSSM parameters. In the 
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where the yield variable is defined as 1^ = n^j/s with being the number density of particle 
X and s the entropy density of the universe. 

Following the procedure given in p!0], we calculate the light-element abundances taking 
account of the hadro-dissociation, photo-dissociation, and p n conversion processes. The 
energy distribution of the final-state particles are calculated by using the HELAS package 
[llj . and the hadronization processes of colored particles are studied by using the PYTHIA 
package [12]. In the Bino-NLSP case, high energy neutrino emitted by the Bino decay may 
scatter off background neutrino and generate energetic e^, which becomes the source of 
energetic photon [13j. In our analysis, we have taken into account the effects of the photo- 
dissociation process induced by photon from the neutrino injection. (However, we found 
that the neutrino-induced processes are less important compared to other processes.) Once 
theoretical values of the primordial light-element abundances are obtained as functions of 
the mass and the lifetime of the NLSP, we compare them with the observed values of the 
primordial abundances. In deriving the constraints on the model, we adopt the following 
observational constraints: 

• D to H ratio [HIIS]: 

(nD/nH)p = (2.82 ± 0.26) x IQ-^ (13) 



^He mass fraction [161 IE] : 



= 0.2516 ± 0.0040. 



(14) 



^He to D ratio [TSlfTOl: 



(n3He/r^D)p < 0.83 + 0.27. 



(15) 



^Li to ^Li ratio ^ [20] 
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(ra6Li/n7Li)p < 0.046 + 0.022 + 0.106. 



(16) 



• ^Li to H ratio [22l[20]: 

logio(ri7Li/nH)p = -9.90 ± 0.09 + 0.35. (17) 
so-called bulk region, the abundance is larger, and is approximately given by 

B VlOOGeV/ 

We have checked that the BBN constraints in such a case are almost the same as the focus-point case. If 
we adopt the abundance in the bulk region, however, Jlp^*^ -' becomes larger than the present dark matter 
density if m^^^ = 100 GeV. Thus we will not consider such a case in the following discussion. 

'^^Asplund et al. reported neLj/nTLj = 0.046 ± 0.022. However, their positive detection has not been fully 
confirmed yet as pointed out in [3T]. Therefore, we consider the observed value as an upper bound. 
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(Here and hereafter the subscript "p" denotes the primordial value inferred by observation.) 
As shown in (fT6|) and (fT7j) . we add positive systematic errors of +0.106 and +0.35 to the 
observational face- values of (n6Li/n7Li)p and logio(^7Li/?T,H)p, respectively. We expect that 
these systematic errors result from possible depletion in stars through rotational mixing ^23j 
or diffusion [211. Since both ^Li and ^Li are destroyed by depletion process, their systematic 
errors are correlated (for more details, see [2Q]). We note here that the standard BBN is 
excluded at more than 4-0" level if we do not adopt the systematic error on ^Li abundance |25j 
(so-called ''Li problem). Thus, to derive a conservative constraint, we add these systematic 
errors. At the end of this letter, we will comment on implications of the u^-DM scenario on 
the ''Li problem!^ 

In Fig. [H we show the constraint from BBN for the Bino-NLSP case on vs. plane 
(with being the lifetime of Bino). The lifetime is related to the fundamental parameters 
via Eq. ([7]); in particular, r^j is proportional to a~^. Taking m^^ = 1.2m^, we calculate a^- 
parameter. In the figure, un-shaded, lightly shaded, and darkly shaded regions indicate the 
region with < 1, 1 < < 10, and > 10, respectively. One can see that the region with 
< 200 GeV is always allowed. This is because, in such a region, the dominant hadronic 
decay processes are four-body ones {B Vfivqq and I'Rlqq'), for which the branching ratio 
is significantly suppressed by the phase-space factor. On the other hand, when the decay 
processes B uruZ and vrIW are kinematically allowed, those three-body decay processes 
have sizable branching ratio, resulting in an enhanced production of hadrons. We can see 
that the lifetime of Bino is constrained to be smaller than < 10^ sec in such a parameter 
region in order not to overproduce deuterium via the hadro-dissociation of ^He. 

In the figure, we also plot contours of constant density parameters. Once the primordial 
abundance of the NLSP is fixed, ' is calculated by using Eq. ([6]). We show the contour 

of f^^'j^'^ '* = Vtc = 0.228 PP; the right-hand side of the line is excluded by the overclosure 
constraint if we adopt the abundance given in Eq. f[T^ . In studying the VR-Dyi scenario, 
we should also consider vr from the MSSM particles in the thermal bath. Following [2], we 
calculate the sneutrino abundance by solving the Boltzmann equation taking account of all 
the relevant sneutrino production processes. The contour of q^,'^'^'^™''') = jg shown in Fig. 
dl the Ojy-parameter is determined by using Eq. ([7]), while the MSSM parameters are taken 
to be rriy^ = 1.2m ^, and hh = 2mg (with hh being the SUSY invariant Higgs mass), ur is 
overproduced below the line of f^^,'^'^''™'*') = with the present choice of parameters. One 
can see that the line is well below the constrained region by BBN. In the present choice 
of parameters, a relatively large value of a^, is needed unless the masses of B and vr are 
degenerate in order to realize j^^,'^'^'^''™^') = However, notice that the relic abundance of 

Ur depends on various parameters. In particular, oj-^^^^™^^) becomes larger when the mass 
difference between vr and vl becomes smaller because the left-right mixing is enhanced. In 

#5 As we will discuss later in considering the ^Li problem, one may adopt a slightly higher value of D to 
H ratio, (nD/nH)p = (3.98lJ]:^?) x lO^^ [14 , and/or that of ^Li to H ratio, logio(n7Li/nH)p = -9.63 ± 0.06 
|26j . We have checked that, even with these observational constraints, the constraints given in Figs. [T] — [3] 
are almost unchanged (as far as the systematic error in the ^Li abundance is taken into account). 
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Figure 1: BBN constraints on the Bino-NLSP case are shown on vs. plane. The un-shaded, 
lightly shaded, and darkly shaded regions are for < 1, 1 < < 10, and > 10, respectively. 
In addition, the contours of 17 = and j^CThermai) _ ^^^^ shown (dotted lines). For 

the calculation of f^?^'^'^™'^'^ we take m^^ = 100 GeV, m£;^ = 1.2m^, hh = 2m^, rrih = 115 GeV, 
and tan /? = 30. 



addition, Q,-^ is also enhanced if we use a larger value of the neutrino Yukawa coupling 
constant; it may happen when we adopt the degenerate neutrino masses. Thus, with other 
choices of parameters, the required value of to realize q^J^^^"^^^^^ = changes. 

Another candidate for the NLSP is the left-handed sneutrino. Such a scenario is attractive 
in the z/^j-DM scenario because the abundance is enhanced if the masses of upt and 

becomes closer. When z/^ is the NLSP, its dominant decay process is ^rZ^*"^ 
and vl —>■ i)jih^*\ Thus, colored and/or charged particles are effectively produced via the 
dominant decay modes. Again, we calculate the light-element abundances taking account 
of the hadro-dissociation, photo-dissociation, and p n conversion processes, and compare 
the resultant light-element abundances with observational constraints given in (fT3|) — (fT7|) . 
The relic abundance of left-handed sneutrino is approximated as [27] : 

^ 2 X 10-1^ X ( \ . (18) 

V 100 GeV/ ^ ' 

The BBN constraints are shown in Fig. [21 We can see that the parameter space is con- 
strained as < 10^ sec (with r,>^ being the lifetime of pl) by the deuterium overproduction 
irrespective of m,>^. This is due to the fact that, if i>l is the NLSP, production of hadrons 
occurs in the dominant decay processes. This is a large contrast to the Bino-NLSP case. 

If pl is the NLSP, its primordial abundance is so small that ^'^^f''^ < as far as 
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Figure 2: BBN constraints on the P^^-NLSP case are shown on rrij;^ vs. r^^^ plane. The un-shaded, 
lightly shaded, and darkly shaded regions are for < 1, 1 < < 10, and > 10, respectively. 
In addition, the contour of qC^^*^™^') = shown in dotted line. Here, we take m£;^ = 100 GeV, 
= 1.2m/>^, fiH = 2m^, ruh = 115 GeV, and tan/5 = 30. 



mo 



< 10 TeV (for 



100 GeV). On the contrary, q^J^^"^"^^^) qq^^ Iqq large as Qc', in 
the figure, we plot the contour of q^^J^'^^^^^^ = Here, the ai,-parameter is determined for 
given values of m£>^ and r;>^, while the MSSM parameters are taken to be = 1.2m/>^, and 
fiH = 2m^. One can see that, when mi>^ > 160 GeV, q^^^'^^'^^^^ = c^n be realized with 
< 10 (which is marginally consistent with the naive order-of-estimate of the a,y-parameter 
in gravity-mediated SUSY breaking scenario). Notice that, even with Oj, ~ 1 (or smaller), 
Q^T^"^"^^^^^ can be large enough if a larger value of is adopted or if i>/j is produced by the 
decay of some exotic particles. 

Next, we consider the f-NLSP case. In this case, because the NLSP is charged, it may 
form a bound state with "^He during the BBN epoch and change the reaction rate [28j. (Such 
an effect is called f-catalyzed effect.) Consequently, ^Li abundance may be significantly 
enhanced if the lifetime of r is longer than ~ 10^ sec. Here, the light-element abundances 
are calculated by including the f-catalyzed effect. Assuming that f is almost right-handed, 
we approximate the primordial abundance as p?] : 



Vf ~ 7 X 10 



-14 



X 



100 GeV J 



(19) 



and calculate the light-element abundance. The numerical result is shown in Fig. [31 As 
in the z/j^-NLSP case, the parameter space Tf > 10^ sec (with being the lifetime of f) 
is excluded. In addition, the ^He is overproduced due to p ^ n conversion process when 
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Figure 3: BBN constraints on tlie f-NLSP case are shown on ruf vs. Tf plane. Tlie un-sliaded, 
lightly shaded, and heavily shaded regions are for < 1^ 1 < < 10, and > 10, respectively. 



In addition, the contour of (^(^T^*^™*^') = ^ shown in dotted line. Here, we take 



m 



100 GeV, 



\.2mf, fifj = 2m^, rrih = 115 GeV, tan/? = 30, and sin^f = 0.3. 



nif > 500 GeV and ~ 10 sec. The ^He constraint becomes more stringent than the 
i>i-NLSP case because the yield variable used in the f-NLSP case is larger. We also show 
the line which satisfies f^?^^''™'^^'' = Vtc, taking = 1.2mf, mj>^ = 1.2mf , hh = 2m^, and 
sin6'f = 0.3. As one can see, the lifetime becomes longer for a given value of compared to 
the case of z/^-NLSF; this is because we have taken a small value of 9r- Even in this case, we 

can see that f](Thermal) _ 

can be realized with < 10 in the parameter region consistent 
with all the BBN constraints. 

Finally we comment on the implication of the z/^j-LSF scenario on the so-called ^Li 
problem. As we have mentioned, the standard BBN is excluded at more than 4-0" level if 
we take the face value of the observational constraints on the ''Li abundance; the theoretical 
prediction of the ^Li abundance becomes significantly larger than the observed value. Even 
though the ^Li problem does not exist if a significant depletion of ''Li occurs in stars, the 
degree of the depletion has not yet been accurately understood. If one adopts models with 
small depletion, the astrophysical or particle-physics solution to the ^Li problem is required. 
It is notable that the ''Li abundance can be reduced if a long-lived particle decays into 
hadrons during the BBN epoch [221 SO]- Thus, in the present case, the decay of the NLSF 
during the BBN may be a solution to the ^Li problem. In the following, we will see that 
the ^Li problem may be solved if B is the NLSF. (For the cases of vl- and f-NLSF, the 
^Li problem is hardly solved because the parameter region with the lifetime longer than 
~ 10^ sec is (almost) excluded, as shown in Figs. [2] and O) 
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To study the Li problem in the present framework, we neglect the systematic error 
(i.e., +0.35 dex) in the observational constraint on ''Li abundance. In addition, because the 
allowed parameter region is sensitive to the observational constraint on '^Li, we consider two 
different observational constraints on ^Li abundance: 

Low ^Li : loglo(^^7Li/^^H)p = -9.90 ± 0.09 [22], (20) 
High^Li: logio(mLi/r^H)p = -9.63 ± 0.06 [26]. (21) 

Notice that the low value corresponds to the one given in (|T7|) . while the high value is from 
measurement using different method to estimate temperature of the atmosphere in dwarf 
halo stars. In addition, because the systematic error in the ^Li to ^Li ratio is correlated to 
that of ''Li, we also remove the systematic error from ( |T6l) : 

(n6Li/n7Li)p < 0.046 + 0.022. (22) 

BBN constraints on the Bino-NLSP case are shown in Fig. |H using the constraints fl20l) 
and ( !22l) (upper panel) or ( !2T|) and ( !22|) (lower panel). Here, constraints on (nD/nH)p, Yp, 
and (n3He/'n-D)p are unchanged from the previous cases; the D to H ratio given in ( TT3i) is 
called "Low D" because of the reason below. As one can see, if we adopt the high value of 
the ^Li to H ratio, all the light-element abundances can be consistent with the observational 
constraints if 10^ sec ^ < 10^ sec. On the contrary, with the low value of ^Li abundance, 
the constraint on the D to H ratio ( fT3l) makes it difficult to solve the ''Li problem. However, 
this conclusion changes if we adopt a slight systematic error in the ^Li abundance, or if a 
different observational constraint on the D to H ratio is adopted. Indeed, in some literature, 
a higher value of the D to H ratio (which is the highest value among the data points for six 
most precise observations [13]) is adopted because D is the most fragile light element and the 
observed values might reflect the abundance after suffering from some destruction processes: 

HighD: (nD/nH)p = (3.981°^?) x 10"^ (23) 

(We call this as "High D.") In Fig. |H we also present the parameter region consistent with 
the constraint fl23l) using the dotted line. As one can see, with fl23|) . the ^Li problem can be 
solved even with the low value of the ^Li to H ratio. 

Notice that, in the parameter region where all the light-element abundances become 
consistent, ^2*7'^'=™='') becomes much smaller than Qc if the constraint ( 120|) or ( 12T1) is adopted. 
However, this fact does not imply that the ur-LSP scenario cannot solve the ^Li problem. 
One possibility is to consider the effects of the decay products of MSSM-LSP after freeze-out; 
indeed, as shown in the figure, f^^^^^ ' — is realized when rrij^ 400 GeV and ~ 10^ sec 
with solving the ^Li problem. 
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